The insulating phases and superfluid-insulator transition of disordered boson chains 
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Using a strong disorder real-space renormalization group (RG), we study the phase diagram 
of a fully disordered chain of interacting bosons. Since this approach does not suffer from run- 
away flows, it allows a direct study of the insulating phases, which are not accessible in a weak 
disorder perturbative treatment. We find that the universal properties of the insulating phase are 
determined by the details and symmetries of the onsite chemical-potential disorder. Three insulating 
phases are possible: (i) an incompressible Mott glass with a finite superfluid susceptibility, (ii) a 
random-singlet glass with diverging compressibility and superfluid susceptibility, (iii) a Bose glass 
with a finite compressibility but diverging superfluid susceptibility. In addition to characterizing 
the insulating phases, we show that the superfluid-insulator transition is always of the Kosterlitz- 
Thouless universality class. 
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Bose systems can be driven into an insulating phase 
by quantum fluctuations due to strong repulsive inter- 
actions and lattice effects. The impact of a disordered 
potential on this superfluid-insulator transition and on 
the nature of the insulating phases is a long standing 
question [1]. In weakly disordered one dimensional sys- 
tems, the momentum-shell renormalization group (RG) 
afforded much progress [2, 3]. But recently, an analy- 
sis using a real space RG suggested that strong disor- 
der can have very different effects on a one dimensional 
Bose system [4]. In particular, it was found that a cer- 
tain type of disorder, that is perturbatively irrelevant 
at weak disorder, can actually induce a transition when 
sufficiently strong, and lead to a new kind of an insula- 
tor termed the Mott glass [4, 5]. The existence of this 
phase transition and the Mott glass were confirmed nu- 
merically [6, 7]. The disorder considered in Refs. [4, 6], 
however, had a very special particle-hole symmetry prop- 
erties that would not be easy to realize in actual experi- 
ments (e.g., [8, 9]). 

In this paper we extend the real-space RG of Ref. [4] 
to treat strong and general disorder potentials, not con- 
fined to the comensurability requirement in Ref. [4] . Our 
starting point is the disordered quantum-rotor model: 



^2 Jj COS (<p j+ l - tfj) , (1) 



This Hamiltonian describes a chain of superfluid grains 
that are connected by a random Josephson coupling Jj 
(see Fig. 1(a)). Each grain has a random charging en- 
ergy Uj , and offset charge Uj , which represents an excess 
screening charge on the site or in its environment. The 
offset charge parameterizes a random on-site chemical po- 
tential [ij — Ufflj. We point out that the lattice model 
(1) can also be derived as a coarse-grained description of 
continuum bosons, where the grain size or lattice spacing 
is set by the healing length of the condensate. 

Using a real-space RG, we show that the system can 
undergo a transition from a superfluid to three possible 



insulating phases, whose nature depends on the symme- 
try properties of the distribution of offset charges, Hj. 
We characterize the insulating phases using the charging 
gap A, the compressibility K = fr and the superfluid 
susceptibility Xs- The latter corresponds to the linear 
response of the order parameter (e 11 ^ ) to the coupling 
jr J2i cos <fi, where the angular brackets and overline de- 
note a quantum expectation value and disorder average 
in that order. The three insulating phases we find are 
illustrated schematically in Fig. 1. They include: (i) 
an incompressible Mott Glass arising for the case of zero 
offset charges Uj = 0, (ii) a glass phase with a diverging 
compressibility which arises if Wj can only take the val- 
ues of or 1/2, and which we term a Random Singlet 
Glass, and (iii) a Bose Glass phase characterized by a 
finite compressibility and a diverging superfluid suscep- 
tibility in the case of a generic distribution of the off- 
set charges, i.e., a non-singular distribution in the range 
— 1/2 < rij < 1/2. We argue that the superfluid phase 
and the nature of the transition are insensitive to the 
disorder properties of the offset charges. 

Our real-space RG analysis, as in random spin chains 
where it was first applied [10-12], eliminates the high- 
est energy scale in the system at each stage through a 
local decimation step. The Hamiltonian then maintains 
its form (Eq. 1), but with renormalized distributions of 
Ji , Ui and Hi . From the flow of the coupling distribu- 
tions we can learn about both the phases of the system 
and its critical points. In particular, the RG provides 
quantitative predictions for the properties of the insulat- 
ing phases and the transition into them, on which we will 
concentrate. 

We proceed to construct the generalized decimation 
procedure. Let us define the global energy scale f2 = 
maxj(Aj, Jj), where Aj = Uj(l — 2\rij\) is the charging 
energy of the site j. For the Hamiltonian in Eq. (1) 
three types of decimation steps are possible. Type 1: 
site decimation. If CI — Aj for some j, we freeze the 
charge on the site j, thus eliminating this degree of free- 
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FIG. 1: The three insulating phases that emerge for different 
classes of disorder, (a) The Mott glass is realized when only 
the Josephson couplings and charging energies are disordered, 
and with no offset charge, n = 0. At large scales it consists of 
effectively disconnected superfluid clusters of random size, (b) 
The "Random singlet" glass appears when the random offset 
charge is restricted to n — 0, 1/2 in terms of the basic boson 
charge. In this phase bosons are delocalized on random pairs 
of remote clusters at all scales, (c) The Bose glass is realized 
for a generic offset-charge distribution. It consists of large 
superfluid clusters acting effectively as weakly coupled spin- 
1/2's in a uniformly-distributed random z-field, given by the 
onsite gap times the sign of n: Asgn(n) — U(l — 2\n\)sgn(n). 



dom. A Josephson coupling Jj-i ~ JjJj+i/Q between 
the sites j — 1 and j + 1 is generated by a virtual tun- 
neling process through the eliminated site (see Ref. [4]). 
Type 2: Bond decimation. If f2 = Jj for some j, sites 
j and j + 1 merge into a superfluid cluster with an ef- 
fective interaction parameter: l/Uj = l/Uj + 1/Uj+i 
(this corresponds to additivity of two capacitances con- 
nected in parallel). The offset charges of the two sites 
simply add up Uj = Uj + TTj+i. Type 3: Doublet for- 
mation. A special RG step is introduced for sites with 
the offset charge rij = 1/2. In this case, if Uj > f2 the 
site j is frozen to its two lowest-lying degenerate charge 
states. Then we can set Uj — > oo and proceed to treat 
the site as a spin- 1/2 degree of freedom (spin-site) with 

-4 sf. The ex- 



p- 



and, similarly, exp(±i(, 



istence of spin sites requires revisiting rules for the bond 
decimation step ( Type 2) . If the strong bond connects 
a spin site to a regular site, then the spin site can be 
simply treated as a half integer site with infinite U. But 
when two spin-sites are strongly coupled by large J cor- 
responding to xy-coupling: ~Jj( s ] s 'j+i + s j s j+i)' they 
freeze into the triplet state with the total spin s = 1 and 
the projection m z = 0: ^= (| Tjij+i) + I tilj+i))- This 
represents a boson resonating in a symmetric superpo- 
sition between the two sites, with a "charging" gap Jj 
to the excited (anti-symmetric) state. Quantum fluctua- 
tions now allow tunneling between sites j — 1 and j + 2 
with strength Jj-iJj+\/ Jj. The last step is identical to 
the real-space RG in the random xx spin chain [12]. 

For a given system, the RG flow, parameterized by 
probability distributions of the couplings, can lead to ei- 



ther a superfluid or an insulating phase. In the super- 
fluid the system coalesces to one large superfluid cluster, 
while in the insulator, it breaks down to clusters with 
large effective charging gaps connected by weak tunnel- 
ing. Quantitative analysis of the RG flow requires in 
general the solution of integro-differential equations for 
the disorder distributions [12]. Remarkably, for the disor- 
der type we consider, the distributions of U, Ji, and Hi, 
are universal in a large vicinity around the fixed points 
of the RG, which govern the superfluid-insulator transi- 
tions. This greatly simplifies the flow equations, and in 
several cases allows an analytical solution. Below we dis- 
cuss the three classes of disorder corresponding to differ- 
ent symmetry properties of the offset-charge distribution. 

No offset charge. The case of rij = 0, for which the 
Hamiltonian (1) is particle-hole symmetric, was analyzed 
in Ref. [4]. For completeness, we review the main results. 
The RG decimation steps (in this case only the decima- 
tions of Types 1 and 2 are needed) translate into the flow 
of the coupling distributions F(U/U) and G(J/Q). Near 
the superfluid-insulator fixed point they acquire the uni- 
versal form: 



Fix) 1=3 exp 



/o(fi) 



, G(x) « goWx"™- 1 . (2) 



Here A is a normalization constant, and x < 1 [13]. Note 
that the typical strength of the Josephson coupling and 
site charging energy monotonically depend on the pa- 
rameters go and fo respectively. In [4] we derived the 
following flow equations for these parameters: 



dfo 
dT 



/o(l - So)i 



dgo 
dT 



-fo 



(3) 



where T = ln(firj/^) and Qq is the initial cutoff energy 
scale. The solutions to these equations are parametrized 
by a single constant C: fo w C + (1 — go) 2 /2. Negative 
C corresponds to the superfluid state in which the charg- 
ing energy is irrelevant (/o flows to zero) and go flows to 
constant larger than unity. Positive C describes the insu- 
lator in which fo is relevant (fo — > oo) and go, indicative 
of the typical Josephson coupling strength, flows to zero. 
The value C — corresponds to the critical point sepa- 
rating the two phases. At this point go flows to 1 and fo 
flows to zero. We note that substituting fo — > y\ makes 
Eqs. (3) assume a standard Kosterlitz-Thouless form. 

The fixed point in the superfluid phase corresponds to 
a classical model with £7j = and a power-law distribu- 
tion of the Josephson couplings with an exponent larger 
than one. The fact that the fixed point is noninteract- 
ing, implies neither the vanishing of the compressibility, 
nor the formation of true long range order: since our 
analysis relies on the grand canonical ensemble, the low- 
est excitation in the superfluid phase corresponds to an 
addition of a particle and not to a phase twist or to a 
Bogoliubov excitation. Thus the vanishing of Uj only 
implies that the energy for adding an extra particle van- 
ishes with the inverse system size, as expected in the 
superfluid phase. Calculation of the compressibility or 
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stiffness of the superfluid requires a more detailed analy- 
sis, that keeps track of the internal Josephson couplings 
and phonon modes within renormalized clusters. 

In this paper we concentrate on the properties of the 
insulating phases, which are particularly interesting since 
they are most drastically affected by the type of the dis- 
order. In the insulating phases, the canonical and grand 
canonical pictures of the excitations are identical, and 
the real space RG correctly describes their properties. 
The insulating phases are best described by a chain con- 
sisting of nearly disconnected clusters, each with its own 
charging gap. The lowest gap corresponds to the energy 
scale fi at which the last site is decimated. In Ref. [4] 
we showed that for the commensurate case with no offset 
charges this gap vanishes with system size L as A <~ 
and is governed by rare and anomalously large superfluid 
clusters. The compressibility vanishes as (In L)/L. In ad- 
dition, the insulator is characterized by a finite superfluid 
susceptibility. We termed this gapless and incompress- 
ible phase a Mott-glass. This phase was also discussed in 
Refs. [5, 14], and confirmed numerically in Refs. [6, 7]. 

Mixed offset: n = (0, 1/2). Let us now allow sites 
to have either zero or half-integer offset charge maintain- 
ing particle-hole symmetry in the problem. Thus the 
Hamiltonian (1) is invariant under the transformation 
rij — > 1 — rij. Such a restriction naturally arises in an 
array of small superconducting grains with a pairing gap 
much larger than the charging energy. The random par- 
ity of the electron number in each grain would lead to 
an offset charge which is randomly either integer or half- 
integer in units of the cooper-pair charge. 

In addition to distributions of U and J, one must now 
follow n's distribution as well, which can be parametrized 
by the three probabilities - q, p, and s - correspond- 
ing to relative densities of integer (rij = 0), half inte- 
ger (rij = 1/2), and doublet (spin-1/2) spin-sites respec- 
tively. On first glance, the problem becomes much more 
complicated since one has to also consider the flow equa- 
tions for p and q (note that p + q + s = 1 due to nor- 
malization). Postponing a detailed description of these 
equations, let us here observe that whether the cluster 
in a renormalized chain is integer or half integer depends 
only on the parity of the number of the bare sites with 
n = 1/2 contained in it. This implies that no matter 
what the fraction of the half-integer sites in the original 
chain was (as long as it is non-zero) large clusters have 
odd or even parity with equal probabilities. For this rea- 
son the distribution quickly flows to a fixed line with 
p = q = (1 — s)/2, provided p > initially. 

This observation significantly simplifies the flow anal- 
ysis, and reduces the number of flow equations to three. 
A detailed analysis of the integro-differential equations 
shows that the distributions of U and J flow to the same 
universal distributions as those described by Eqs. (2), 



and we obtain a simplified set of flow equations: 
^ = / [i- ff0 (i- *)(! + /„)], 
^ = -f [(l- S )/o + 2^ ], (4) 

ar = y(i-0 -»-(!--)■ 

These have two different families of solutions, corre- 
sponding to two different phases. In one family, the flow 
is to a stable fixed line with s — 0, / = and g > 1: a 
superfluid phase identical to that of the no offset charge 
case. Close to this line, the flow equations reduce to those 
for the integer case Eqs. (3), except for an extra factor of 
1/2 in the equation for go. This factor enters because only 
integer sites with a charge gap can renormalize J. The 
fact that the random offset is not important for describ- 
ing the superfluid is not surprising: large local particle 
number fluctuations in this phase completely overwhelm 
small fluctuations in the offset charge. 

The superfluid-insulator transition is also very simi- 
lar to the n = case, and it is described by essentially 
the same Kosterlitz-Thouless critical point with go = 1, 
/o = and s = 0, similar to Eqs. (3). In the insulating 
phase at go < 1, however, f becomes relevant, as well as 
the spin-site density s, which quickly flows to s = 1. This 
second family of solutions of (4) describes a completely 
different insulating phase than in the case of zero offset 
charges; it corresponds to an effective spin-1/2 chain with 
random x — y ferromagnetic couplings, a model that was 
analyzed in detail in Ref. [12]. Its ground state consists of 
random non-crossing pairs of sites at varying distances, in 
which the spins form the m z = state, (|f|) + lit))- 
This phase is termed the Random- Singlet Glass, reafirm- 
ing the connection with anti-ferromagnetic random spin 
chains [15, 16]. In the bosonic language the ground state 
has bosons delocalized randomly between pairs of sites. 

Many properties of the Random-Singlet Glass can be 
found directly from the analysis of Ref. [12]. The en- 
ergy scale associated with breaking a singlet between 
sites of distance I is fig = flocxp(—\/£). By setting 
I = L, the system size, we obtain the scaling of the 
gap vs. the system size. Following the identification 
rii = i + sf, the compressibility k and superfluid sus- 
ceptibility Xs m the insulating phase correspond to the 
susceptibilities of a random spin chain to a Zeeman field 
applied in the z and x directions respectively. From Ref. 
[12] we see that both susceptibilities diverge at the limit 
of small n as - Xs (ty ~ l/fllog 3 (fl /fl). While 

the superfluid stiffness vanishes in the thermodynamic 
limit, unlike the stiffness of the Mott glass, it vanishes 
only sub-exponentially with system size: p s oc e~^. 
In addition, note that go, indicative of the strength of 
the Josephson coupling, flows to zero as go ~ 1/r (as 
results from Eq. (4) with s = 1), much slower than 
go ~ exp[e~ r ] = exp[— 1/fi] in the Mott glass. 

Generic chemical potential disorder. When all 
possible offsets — 1/2 < rij < 1/2 are allowed, the rele- 
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vant energy scale for the RG site decimation is not U, but 
rather the local gap Aj = Ui(l — 2|nj|). The interaction 
Ui is allowed to exceed il by an amount which depends 
on the local charge offset: 



this disorder average from the distribution of Aj: 



Uj < n/(l - 2\Hj 



(5) 



Thus we must consider a joint distribution for U and n, 
which makes the RG flow equations rather cumbersome. 
Their analysis, however, reveals a rather intuitive behav- 
ior. Below we describe the key steps of the derivation. 
The first step is to note that n disorder width is a 



relevant variable due to the rule 



ij+i 



under 



the decimation of a strong bond. This implies that as 
the effective sites grow with the RG, their offset charges 
quickly become uniformly distributed between —1/2 and 
1/2, i.e., the largest disorder allowed. This observation 
significantly simplifies the analysis. In particular, it is 
straightforward to check that the distributions of U and 
J again approach the universal functions (2). In the su- 
pcrfluid regime and at the transition point the flow is 
also governed by the equations (3) with an extra factor 
of one half in the equation for g - Thus the system with 
generic disorder undergoes the same transition as in the 
two cases discussed above. Similarly, the elementary un- 
derstanding of this result is that for small interactions 
the bounded disorder in the offset charge is overwhelmed 
by large particle number fluctuations. In the insulating 
side, /o flows to large values; the joint distribution of U/fl 
and n is still given by (2): F joint (U/fl,n) = F(U/fl), 
but with the upper bound of U/Q obeying Eq. (5). At 
large / , this form of F jomt (U/Q,n) leads to the uni- 
form distribution of the charging gaps A of each site: 
-ff(A) w 1/17. Most importantly, H(A) is non-singular 
at A -> 0. 

Deep in the insulating phase, since the typical charging 
energy U is large, each site can be treated as a doublet of 
states which is split locally by the energy A < 0, while 
all other states lie at energies above the cutoff. When 
A = this doublet is exactly degenerate and represents 
the spin-1/2 degree of freedom discussed above. Thus 
the spin-1/2 description is also applicable to the case of 
nonzero A with the latter playing the role of the Zeeman 
field along the z-axis. The effective Hamiltonian describ- 
ing the chain becomes 



X X 

S 3 S 



S 1 jS V 



3 + 1 J 



3 



(6) 



where sf = rij — 1/2, and here Aj = Uj(l — 2\nj\)sgn(rij) 
can be both positive and negative (its distribution is 
H(Asgn(nj) = 1/20). Because Jj <!lwe can calcu- 
late many properties of this phase by considering a single 
site subject to a random Zeeman field. First, as implied 
by the gap distribution, in this phase the energy-length 
scaling is Ql ~ 1/L = P- Next, the compressibility is 
given by the resp onse to an external field in the z direc- 
tion: k = pd(s z ) /dh z ext — 2H(0) = p/Q = k , a constant 
at low energies. The superfluid sus cepti bility is the re- 
sponse to a transverse field: \s = pd(s x ) /dh x xt . We find 



Xs 



dA 



log 



n 2 ^ ViK.t) 2 + A 2 2 »\\h 



(7) 



which diverges as h x xt — > 0, with a functional form result- 
ing from the non-singular behavior of H(A). In a finite 
chain this divergence is cut off by the smallest |Aj|, and 
leads to \s ~ log£. These properties, namely finite com- 
pressibility and diverging superfluid susceptibility, coin- 
cide with the Bose Glass phase discussed in Rcf. [1]. 



n disorder 


Glass type 


A 


K 


Xs 




0, 1/2 


Mott 
Random- 
singlet 


i 

logL 

e -cvT 


log L n 
L u 

1 cVT 


const 

1 CsfL 


-1/2 <n< 1/2 


Bose 


1 
L 


K 


f log-L 



TABLE I: The gap A, compressiblity k, and SF susceptibility 
Xs, of the insulating phases realized for the different classes 
of disorder in the offset charge n. Here L is the system size 
and c denotes a nonuniversal constant. 



In summary, we extended the real-space RG analysis of 
Rcf. [4] to address non-commensurate random chemical 
potential. We found that the symmetry of the Hamil- 
tonian and the details of the disorder, as encoded in 
the distribution of the offset Uj, do not affect the uni- 
versal properties of the superfluid phase and superfluid- 
insulator transition, but at the same time it completely 
determines the insulating phase of the system. In addi- 
tion to the Mott-Glass of the integer filling case, we find 
the Random- Singlet Glass for the mixed n = 0, 1/2 case, 
and for generic chemical-potential disorder, a phase we 
identify as the Bose-Glass[l]. We note that if we restrict 
n to a set of rational numbers, we would still obtain ei- 
ther the random-singlet glass, or the Mott-glass. All our 
results are supported by a numerical real-space RG study 
to be published separately. 

Our focus here was the universal thermodynamic prop- 
erties of the insulating phases, as summarized in Table 
I. An outstanding question, which we leave to a future 
publication, is the effective Luttingcr parameter at the 
critical point, and how it compares to the low-disorder 
motivated results of Rcf. [2], predicting K = ^Jnp s = 3/2. 
Our initial results, however, indicate that the critical 
value of the Luttinger parameter at sufficiently strong 
disorder is non- universal, and exceeds 3/2; this possibil- 
ity does not contradict a general thermodynamic state- 
ment of Ref . [6] . 
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